Nonlinear and chaotic resonances in solar activity 
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It is shown that, the wavelet regression detrended fluctuations of the monthly sunspot number 
for 1749-2009 years exhibit strong periodicity with a period approximately equal to 3.7 years. The 
wavelet regression method detrends the data from the approximately 11-years period. Therefore, 
it is suggested that the one-third subharmonic resonance can be considered as a background for 
the 11-years solar cycle. It is also shown that the broad-band part of the wavelet regression de- 
trended fluctuations spectrum exhibits an exponential decay that, together with the positive largest 
Lyapunov exponent, are the hallmarks of chaos. Using a complex-time analytic approach the rate 
of the exponential decay of the broad-band part of the spectrum has been theoretically related to 
the Carrington solar rotation period. Relation of the driving period of the subharmonic resonance 
(3.7-years) to the active longitude flip-flop phenomenon, in which the dominant part of the sunspot 
activity changes the longitude every 3.7 years on average, has been briefly discussed. 

PACS numbers: 05.45.a, 47.65. Md, 96.60.qd 



in 

I • 
O 



(N 
> 

m 

o 
o 



X 



The solar activity is chaotic but has a well-defined 
mean period of about 11 years. The 11-year cycle is well 
known for more than a century and a half. Despite this, 
nature of the 11-year cycle is still unknown. Most of the 
regression methods are linear in responses and statisti- 
cal analyses of the experimental sunspot data was domi- 
nated by linear stochastic methods, while it was recently 
rigorously shown in Ref. [l[ that a nonlinear dynamical 
mechanism (presumably a driven nonlinear oscillator) de- 
termines the sunspot cycle. Figure 1 shows the monthly 
sunspot number (dashed line) for the period 1749-2009 
years (the data were taken from Ref. 0). The solid 
curve (trend) corresponds to a wavelet (symmlet) regres- 
sion of the data (cf. Ref. Q). Figure 2 shows correspond- 
ing detrended fluctuations, which produce a statistically 
stationary set of data. At the nonlinear nonparametric 
wavelet regression one chooses a relatively small number 
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FIG. 1: The monthly sunspot number (dashed line) for the 
period 1749-2009 years. The data were taken from Ref. 01 . 
The solid curve (trend) corresponds to a wavelet (symmlet) 
regression of the data. 



of wavelet coefficients to represent the underlying regres- 
sion function. A threshold method is used to keep or 
kill the wavelet coefficients. In this case, in particular, 
the Universal (VisuShrink) thresholding rule with a soft 
thresholding function was used. At the wavelet regres- 
sion the demands to smoothness of the function being 
estimated are relaxed considerably in comparison to the 
traditional methods. Figure 3 shows a spectrum of the 
wavelet regression detrended data calculated using the 
maximum entropy method (because it provides an opti- 
mal spectral resolution even for small data sets). In Fig. 
3 one can see a well defined peak corresponding to period 
~ 3.7 years. The wavelet regression method detrends the 
data from the approximately 11-years period (cf. Fig. 1). 
Therefore, it is plausible that the one-third subharmonic 
resonance can be considered as a background for the 
11-years solar cycle: 11/3.7 ~ 3. Indeed, it is known 
Q that interaction of the Alfven waves (generated in a 
highly magnetized plasma by a cavity's moving bound- 
aries) with slow magnetosonic waves can be described 
using Duffing oscillators (see also Refs. [l|,[ll). Let us 
imagine a forced excitable system with a large amount of 
loosely coupled degrees of freedom schematically repre- 
sented by Duffing oscillators (which has become a classic 
model for analysis of nonlinear phenomena and can ex- 
hibit both deterministic and chaotic behavior de- 
pending on the parameters range) with a wide range of 
the natural frequencies loq : 



X + uJqX + jx + /3x'^ ~ F sincji 



(1) 



where x denotes the temporal derivative of a;, /3 is the 
strength of nonlincarity, and F and w are characteristic 
of a driving force. It is known (see for instance Ref. Q) 
that when uj « 3a;o and /3 <^ 1 the equation (1) has a 
resonant solution 
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FIG. 2: The wavelet regression detrended fluctuations from 
the data shown in Fig. 1. 




FIG. 4: The same as in Fig. 3 but in semi-logarithmical scales. 
The dashed straight line indicates an exponential decay. 
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FIG. 3: Spectrum of the wavelet regression detrended fluctu- 
ations shown in Fig. 2. 



where the amphtude a and the phase ip are certain 
constants. This is so-called one-third subharmonic 
resonance with the driving frequency u corresponding 
approximately to 3.7 years period (the peak in Fig. 3 
corresponds to the second term in the right-hand side 
of the Eq. (2) while the first term has been detrended). 
For the considered system of the oscillators an effect of 
synchronization can take place and, as a consequence 
of this synchronization, the characteristic peaks in the 
spectra of partial oscillations coincide [8j. It can be 
useful to note, for the solar activity modeling, that the 
odd-term subharmonic resonance is a consequence of the 
reflection symmetry of the natural nonlinear oscillators 
(invariance to the transformation x —x). 

In order to understand appearance of the 3. 7- years pe- 
riod let us represent the spectrum shown in Fig. 3 in 
semi-logarithmical scales: figure 4. In these scales an ex- 
ponential behavior corresponds to a straight line. It is 



known, that both stochastic and deterministic processes 
can result in the broad-band part of the spectrum, but 
the decay in the spectral power is different for the two 
cases. The exponential decay indicates that the broad- 
band spectrum for these data arises from a deterministic 
rather than a stochastic process. For a wide class of 
deterministic systems a broad-band spectrum with expo- 
nential decay is a generic feature of their chaotic solu- 
tions Refs. Q-fl^l- A wavy exponential decay (see Fig. 
4) is a characteristic of a chaotic behavior generated by 
time-delay differential equations [10]. A classic example 
of time-delay differential equation with chaotic solutions 
is the Mackey-Glass equation: 
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Figure 5 shows spectrum of a solution of this equation for 
the time-delay r = 30. The dashed straight line indicates 
an exponential decay (cf. Fig. 4). 

In the dynamo models that have physically distinct 
source layers the finite time is required in order to trans- 
port magnetic flux from one layer to another [l^, it is 
especially significant for those dynamo models that have 
spatially segregated source regions for the poloidal and 
toroidal magnetic field components (such as, for instance, 
the Babcock- Leighton dynamo mechanism [l3|)- In the 
global dynamo models that include meridional circula- 
tion the time delay related to the circulation should be 
comparable to global rotation period (see below). 

Nature of the exponential decay of the power spectra 
of the chaotic systems is still an unsolved mathematical 
problem. A progress in solution of this problem has been 
achieved by the use of the analytical continuation of the 
equations in the complex domain (see, for instance, 12|). 
In this approach the exponential decay of chaotic spec- 
trum is related to a singularity in the plane of complex 
time, which lies nearest to the real axis (see the insert in 
Fig. 5). Distance between this singularity and the real 
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FIG. 5: Spectrum of the Mackey-Glass chaotic time-series. 
The dashed straight hne indicates an exponential decay. The 
insert shows a sketch of corresponding complex-time plane. 



axis determines the rate of the exponential decay. For 
many interesting cases chaotic solutions are analytic in a 
finite strip around the real time axis. This takes place, 
for instance for attractors bounded in the real domain 
(the Lorentz attractor, for instance). In this case the ra- 
dius of convergence of the Taylor series is also bounded 
(uniformly) at any real time. 

Let us consider, for simplicity, solution u(t) with simple 
poles only, and to define the Fourier transform as follows 

u{f) = (27r)-i/2 / dt e-'2'^^*?i(i) (2) 

J-Te/2 

Then using the theorem of residues 

uif) = z(27r)i/2 ^ exp(z2^/x, - \27r fy,\) (3) 



where Rj are the poles residue and Xj + iyj are their lo- 
cation in the relevant half plane, one obtains asymptotic 
behavior of the spectrum E{f) — |u(/)p at large / 



E{f) ~ exp(-47r|y„,„| /) 



(4) 



where j/mm is the imaginary part of the location of the 
pole which lies nearest to the real axis. In the case of 
symmetric analytic strip with a width A = 2\ymin\'- 



E{f ) ^ exp(-27rA /) 



(5) 



(cf. the insert in Fig. 5). 

The chaotic spectrum provides two different char- 
acteristic time-scales for the chaotic system: a period 
corresponding to fundamental frequency of the system, 
ffund, and a period corresponding to the exponential 
decay rate, 27rA (cf. Eq. (5)). The fundamental period 
can be estimated using position of the low-frequency 
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FIG. 6: The pertaining average maximal Lyapunov exponent 
at the pertaining time, calculated for the same data as those 
used for calculation of the spectrum (Figs. 3 and 4). The 
dashed straight line indicates convergence to a positive value. 



peak (cf. Figs. 4 and 5), while the exponential decay 
rate period 27rA can be estimated using the slope of the 
straight line of the broad-band part of the spectrum 
in the semi-logarithmic representation. In the case of 
the global solar dynamo the width of the analytic strip 
A can be theoretically estimated using the Carrington 
solar rotation period: A ~ Tj, ~ 25.38 days. This period 
roughly corresponds to the solar rotation at a latitude 
of 26 deg, which is consistent with the typical latitude 
of sunspots (cf. Fig. 4). 

Additionally to the exponential spectrum (Fig. 4), let 
us check the chaotic character of the wavelet regression 
detrended fiuctuations calculating the largest Lyapunov 
exponent: Xmax- A strong indicator for the presence of 
chaos in the examined time series is condition Xmax > 0. 
If this is the case, then we have so-called exponential 
instability. Namely, two arbitrary close trajectories of 
the system will diverge apart exponentially, that is the 
hallmark of chaos. To calculate Xmax we used a direct 
algorithm developed by Wolf et al. [ll]. Figure 6 shows 
the pertaining average maximal Lyapunov exponent at 
the pertaining time, calculated for the data set shown in 
Fig. 2. The largest Lyapunov exponent converges very 
well to a positive value Xmax — 0.286 mon~^ > 0. 

It should be noted that the same period ^ 3.7 
years was recently found for the so-called fiip-flop 
phenomenon of the active longitudes in solar activity 
|16| , ITf . Sunspots are tend to pop up preferably in 
certain latitudinal domains and move toward the equator 
due to the 11-year cycle. Recently, strong indications 
of non-uniform longitudinal distribution of sunspots 
(active longitudes) was reported and analyzed in a 
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dynamic frame related to the mean latitude of sunspot 
formation, in whicli the active longitudes persist for the 



last eleven solar 11-years cycles (see Refs. [16[,[17[ and 
references therein). At any given time, one of the two 
active longitudes (approximately 180° apart) exhibits 
a stronger activity - dominance. Observed alternation 
of the active longitudes dominance in 3.7 years on 
average was called as flip-flop phenomenon [l6|. It 
seems rather plausible that the observed flip-flop period 
and the period of the wavelet regression detrended 
fluctuations of solar activity (Fig. 4) have the same 
origin. In this vein, the observation [lB|,[l3 that the 
period of the flip-flop phenomenon follows to variations 
of the real length of the sunspot cycle (which has the 
11-years period on average only) supports the idea of 
the one-third subharmonic resonance as a background of 
the 11-years cycle of solar activity. 

The author is grateful to SIDC-team, World Data Cen- 
ter for the Sunspot Index, Royal Observatory of Belgium 
for sharing their data. A software provided by K. Yosh- 
ioka was used at the computations. 
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